Confined water in the low hydration regime. 
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Abstract 

Molecular dynamics results on water confined in a silica pore in the low hydration regime are presented. Strong layering 
effects are found due to the hydrophilic character of the substrate. The local properties of water are studied as function of both 
temperature and hydration level. The interaction of the thin films of water with the silica atoms induces a strong distortion 
of the hydrogen bond network. The residence time of the water molecules is dependent on the distance from the surface. Its 
behavior shows a transition from a brownian to a non-brownian regime approaching the substrate in agreement with results 
found in studies of water at contact with globular proteins. 

PACS numbers: 61.20.Ja, 61.20.-p, 61.25.-f 



I. INTRODUCTION 

Thin films of water confined in restricted geometries 
play an important role in many different phenomena in 
geology, chemistry, biology. Important processes taking 
place in minerals, plants, biological membranes are de- 
termined by the presence of water confined in narwupares 
inside the system or at contact with the surface .tyflEl In 
biology for instance the enzymatic activity of globular 
proteins or the function of bioiogical membranes are due 
to the interaction with wateitl, and the decrease of the. 
hydration level may inhibit their biological functionality.!!! 
It is still difficult to predict accurately how the properties 
of water are modified when it is at contact with the var- 
ious substrates or confined in different geometries. Nev- 
ertheless recent theoretical and computational work on 
water confined in hydrophobic rigid environments ind 
cate strong changes in its thermodynamical behaviorjj' 

A strong reduction of the tetrahedral order with conse- 
quent distortion of the hydrogen bond network has been 
evidenced by neutron scattering experiments with iso- 
tonic- substitution performed on water confined in Vy- 
coid'EI, by experiments performed on freezing ofpsvater 
confined in an hydrophilic environment of micaa and 
in water confined in vcrmiculite clay.lij These findings 
are iH-,aspeement with recent computer simulation re- 
sults. OII3II3 Particularly relevant are the observations in 
experiments and computer simulation of the formation of 
distinct layer of water with different structural and dy- 
namical properties: an interfacial layer close to the sub- 
strate, where water molecules show a slow dynamics, and 
a layer far ApiB the, surface with almost bulk like prop- 
ertiesMliyOBNea The layering effect is expected to 
be relevant in many phenomena associated with confined 

relax- 
while 



water. There are experimental evidences 
ation of water at contact with proteinsEir 
signatures of a similar behavior has been found in water 
confined in Vycor for low hydrationsHj 

As long as a detailed description of the microscopic 
properties of water at interface with complex materials 
like biological macromolecules is difficult to be obtained 



it could be appropriate to use model systems with well 
defined geometry and interaction between water and sub- 
strate. We present here the results obtained by a molec- 
ular dynamics (MD) study of a model of water confined 
in Vycor at low hydration levels. By computer simula- 
tion we can freely vary temperature and hydration level 
and make an accurate study of the changes induced in 
the properties of water by the confinement. This work is 
part of a detailed study of the. sijpjj^ital and dynamical 
properties of confined water Jij'E^lllSEalHlj 

The low hydration regime is somehow more important 
to be explored at least for two main reasons: (i) the ef- 
fect of the substrate on the confined water is expected to 
be enhanced; (ii) the experimental observations are the 
result of an averaging over different water layers and an 
accurate analysis of the low hydration regime is the only 
way to gain some insight on the water-substrate interac- 
tion. 

An accurate layer analysis enables us to examine the 
modifications of the main structural properties of con- 
fined water and to make connection with recent exper- 
imental and computer simulation results which concern 
water at contact with proteins. 

In Sec. II we briefly describe our simulation method, 
in Sec. Ill, IV and V we present the main results of our 
calculation. Sec. VI is devoted to conclusions. 



II. COMPUTER SIMULATION OF CONFINED 
WATER 

The computer simulation is performed in the micro- 
canonical ensemble on water molecules described by the 
simple point charge/extended (SPC/E) site model poten- 
tial inserted in a cylindrical cavity of radius 20 A. The 
cavity has been carved in a cubic cell of silica glass ob- 
tains 
work.' 1 



„ ..computer simulation, as described in previous 
EjEJ The system is modeled to represent water 
confined in Vycor, the internal surface of the cylindrical 
cell is corrugated and composed by silicon atoms (Si), 
bridging oxygens (BO) which are bonded to two silicons 



FIG. 1: snapshots of the configuration of water molecules at 
room temperature: from above and from the left Nw = 500, 
N w = 1000, N w = 1500. 

and non-bridging oxygens (NBO) connected only to one 
silicon. The NBO are saturated with acidic hydrogens 
(AH) . The interaction of water with the substrate is mod- 
eled by assigning different charges to the sites represent- 
ing the atoms of Vycor. A Lennard- Jones interaction is 
assumed between the oxygens-af water and the BO and 
NBO of the confining system£3 Periodic boundary con- 
ditions are applied along the axis of the cylinder assumed 
as the z direction. The motion is confined in the xy plane, 
where the distance from th e axis is defined in terms of 
the radius R = \J x 2 + y 2 . During the simulation the 
substrate is kept rigid. _ 
The density of the full hydration in the experimentsQ 
0.0297 A~ 3 corresponds in our system to roughly Nw = 
2600 molecules. We performed our simulation at three 
different hydration levels Nw = 500 which corresponds 
to 19 % of hydration, N w = 1000 (38 %) and N w = 1500 
(58 %). In Fig. |l| the snapshots of the configurations at 
the three different hydration levels show the hydrophilic 
character of the substrate. 



III. DENSITY PROFILES 

In the density profiles reported in Fig. ^ it is evi- 
dent for all the hydrations the formation of a well de- 
fined first layer of molecules close to the surface, while 
a second layer grows up with the hydration level. Few 
molecules penetrate and are trapped inside the substrate 
(R > 20 A). We notice that at the highest hydration 
(Nw — 1500) the density in the double layer structure 
increases above the value of density of water at ambient 
conditions. The decrease of temperature has very little 
effect on the density profiles, apart from the sharpening 
of the second layer at the hydration Nw = 1500. 



FIG. 2: Density profiles of the confined water along the pore 
radius for different hydration levels: from below Nw = 500, 
Nw = 1000, Nw = 1500 at two different temperatures: T = 
300 K (solid line), T = 220 K (dotted line). The density 
is normalized to the value for water at ambient temperature 
p bulk = 0.0334 A" 3 . 

It is clear from Fig. [I] that the arrangement of the 
molecules close to the Vycor surface is not uniform and 
the formation of clusters of different sizes is observed, 
this is even more evident in Fig. |^ where we report an 
histogram of the distribution of the molecules in the xy 
plane along a circumference at distance R = 17.5 A where 
the density profile reaches the maximum. The cluster 
sizes fluctuate up to hundred molecules with an average 
value of 40 molecules, this is in agreement with neutron 
diffraction experimental] where clusters of the same size 
are observed at similar value of the hydration. 



IV. LAYER ANALYSIS OF THE HYDROGEN 
BOND 

Figure ^ shows the average number of hydrogen 
bond(HB) per molecule along the pore radius for the dif- 
ferent hydration levels. We assume that two neighbors 
water molecules are luairogen bonded if the H — O ■ ■ ■ O 
angle is less than 30° Bj Two water molecules are consid- 
ered neighbors if their 0-0 separation is less than 3.35 A. 
The number of water- water HB N~j£ (R) decreases ap- 
proaching the surface, while inside the double layer struc- 
ture it grows up the number of HB of the water with 
the Vycor surface N^(R). The water molecules form 
preferentially hydrogen bonds with the bridging oxygens. 
The arrangement of the clusters close to the Vycor sur- 
face is mainly determined by this type of bonds. The 
number of HB between water molecules and the acidic 
H atoms is negligible on the scale of the figure. The 
N£B(R) curves do not show substantial differences as 
function of hydration level of the pore. 

In Fig. H the radial profile of the water- Vycor bonds 
is compared with the density profile of the bridging 
oxygens at the surface; by considering that the bond 
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FIG. 3: Instantaneous density profile for Nw = 500 at room 
temperature in the xy plane along a circumference of radius 
R — 17.5 A corresponding to a maximum in the average den- 
sity profile at the same hydration (see Fig. 2). The density 
profile is obtained from the snapshot in Fig. 1 and shows the 
arrangement of the molecules in clusters. 



length is approximately 2 A there is a clear correspon- 
dence between the two profiles. It is evident that the 
bridging oxygens substitute as acceptors the water oxy- 
gens approaching the surface since the total number of 
HB N% t B (R) = N*™(R) + N£ V B (R) , also reported in 
Fig. |J, remains practically constant. The maximum num- 
ber of water-water HB grows with increasing hydration 
level since more water molecules become available for the 
formation of bonds. The maximum is located around 
R = 16 A where the density profile shows a minimum 
(see Fig. ||), in agreement with the fact that the hydrogen 
bonds increase with decreasing, density at room temper- 
ature and upon supercooling.ta 

At lower temperatures the number of water-water and 
water- Vycor hydrogen bonds increases but the profiles 
are similar to those shown in Fig. [| at T — 300 K. 

In Fig. H we report the radial profiles of the average 
number of water-water hydrogen bonds compared with 
the profiles of the average number of neighbors N nn . The 
ratio between the two quantities N BB (R) /N nn (R) esti- 
mates the number of molecules forming hydrogen bonds 
between the neighbors. It can be considered as a mea- 
sure of the rigidity of the hydrogen bond network and it 
is reported in Fig. [?] for selected values of the hydration 
and temperature. The ratio is smaller close to the surface 
(R > 16 A) where the hydrogen bonds are made mainly 
with the Vycor surface. Below 16 A each molecule is able 
to make hydrogen bonds with almost all its neighbors. 

The network becomes more rigid on lowering the tem- 
perature (see Fig. ^) since the number of hydrogen bonds 
increases while the number of neighbors is almost insen- 
sible to the variation of temperature, like in the density 
profiles in Fig. ||. This might be possibly related to the 
recently proposed fragile to strong transition for strongly 
confined water .113 

We can conclude that the water molecules rearrange 




FIG. 4: Variation of the average number of hydrogen 
bonds per molecule along the pore radius at T — 300 K for 
N w = 500 (top), N w = 1000 (middle) and N w = 1500 (bot- 
tom) at room temperature. The solid line is the water-water 
Nww(R)y the dotted line is the water- Vycor N^ (R) and the 
long dashed line represents the total N^(R). 



themselves at decreasing temperature in order to maxi- 
mize the number of hydrogen bonds. We observe more- 
over that the rigidity increases at decreasing hydration 
(see inset of Fig. [?]) with some important consequences 
for the residence times, as discussed later. 

We investigated how the short range order present 
in bulk liquid water is modified upon confinement by 
looking at the distribution of the hydrogen bond angles. 
Fig. U shows the distribution of cos(9), where 9 is the 
angle between the intramolecular bond O-H and the in- 
termolecular O- • O vector, for Nw — 1000 and different 
temperatures. The water-water linear hydrogen bonds 
are still favored like in bulk water, there is however a 
more pronounced tail toward large angle with respect to 
bulk water. By lowering the temperature the distribu- 
tion peak becomes sharper indicating an increase in the 
rigidity of the bonds. The distortion in the angular dis- 
tribution with respect to the bulk is due to the water 
molecules close to the pore surface, as clearly appears 
in the inset of Fig. B where for Nw = 500 the sepa- 
rated contributions from the layer < R < 16 A and 
16 < R < 20 A are shown. For the molecules close to 
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FIG. 5: The profile of the average number of water- Vy cor hy- 
drogen bonds per molecule along the pore radius (top) com- 
pared with a profile of the number of bridging oxygens (bot- 
tom) on the Vycor substrate (N w = 1000 and T = 300 K). 

the surface the peak of the angular distribution is de- 
pressed while there is a large increase of the tail at high 
angle, the angular distribution of the molecules closer to 
the interior of the pore is similar to the one of the bulk 
water. 

The distortion of the hydrogen bonds is confirmed by 
experimental studies on wateri-confined in VycorEIH and 
at contact with other surfacesfl 

We studied by layer analysis the local arrangement of 
the nearest neighbor molecules by looking to the distri- 
bution of the angle 7 between the two vectors joining 
the oxygen atom of a water molecule with the oxygen 
atoms of the two closest water neighbors. This distribu- 
tion in the computer simulation of bulk water at ambi- 
ent condition has a well defined peak around 7 w 109° 
(cos(7) ~ —0.326), this peak is a signature of the tetrahe- 
dral order present in liquid water at short range distance. 
There is also a secondary peak at 7 w 54° correspond- 
ing to interstitial neighbor molecules. These findings are 
in agreement with, experimental results on water at am- 
bient conditions Ea In Fig. ^| we report a layer analysis 
of the distribution function of cos(7) of water in Vy- 
cor for different hydrations. The molecules in the layer 
< R < 16 A for N w = 1500 and N w = 1000 show 
a distribution very similar to what found in the simula- 
tion of bulk water, the distributions become sharper at 
lower temperature. For Nw — 500 the peak shifts to- 
ward 90°, the value at which the maximum is located for 
T = 220 K. This behavior for the lowest hydration is in 
agreement with experimental results indicating that in 
confined water freezing in a cubic crystalline phase could 



FIG. 6: Comparison of the water-water hydrogen bonds per 
molecule N^ w (R) (top) and the average number of nearest- 
neighbor water molecules N nn (R) (bottom) along the pore 
radius for the hydration level Nw = 1500, Nw = 1000 and 
Nw = 500 in descending order at temperatures T = 300 K 
(solid), T — 220 K (long-dashed). 
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FIG. 7: Ratio N%£ (R)/N nn (R) along the pore radius for 
N w = 1000 at temperatures T — 300 K (solid), T = 220 K 
(long-dashed). In the inset N ww (R) / N n n(R) is reported at 
T = 300 K for the hydrations Nw = 1500 (solid line), N w = 
1000 (dotted line) and Nw = 500 (long dashed line). 



be favored S3 

As a matter of fact a good short range order is present 
for all the hydrations in the layer < R < 16 A where 
the water molecules are able to form hydrogen bonds with 
almost all their neighbors as discussed above (see Fig. ||). 

In the layer closest to the Vycor surface 16<i?<18A 
the tetrahedral arrangement is completely lost due to the 
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FIG. 8: Hydrogen bond angle distribution function, P(cos9) 
(see text for definition) for Nw = 1000 and different temper- 
ature: T = 300 K (solid) ,T = 260 K (dashed), T = 220 K 
(long dashed). In the inset for Nw = 500, layer analysis of 
the distribution at T=300, < R < 16 A (long dashed), 
16 < R < 20 A (solid). 



formation of hydrogen bonds between water molecules 
and the substrate. The temperature has very little effect 
on these distributions. In the layer 18 < R < 20 A the 
specific shape of the distribution is not very significant 
and indicates the strong distortion of the hydrogen bonds 
formed by the molecules found in the inner surface of the 
pore. 
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FIG. 9: Layer analysis of the angular distributions of the angle 
between the oxygens of three nn water molecules, shown for 
N w = 500 (top), Nw = 1000 (middle), N = 1500 (bottom) at 
T = 300 K (solid line) and T = 220 K (dashed line). For each 
hydration the figure shows the angular distribution for the 
layer < R < 16 A on the left, the layer 16 < R < 18 A in the 
center and the layer 18 < R < 20 A closest to the substrate 
on the right. For comparison the angular distribution for 
ambient water is shown as a dotted line in the first inset at 
bottom on the left. 



V. RESIDENCE TIME 

The dynamics of the molecules of confined water is 
strongly determined by the interaction with the sub- 
strate. In Fig. 10 is reported the trajectory of the center 



of mass of a molecule which resides at the beginning in 
the first layer at roughly 2 A from the surface for the 
case Nw = 500. At ambient temperature the molecule 
oscillates mainly between the two layers and persists in 
a position for fractions of nanoseconds, at lower tem- 
perature the molecule oscillates for all the time explored 
around the same position and does not go out from its 
initial layer. 

Fig. [n] shows how the average residence time of the 
water molecules changes along the pore radius for differ- 
ent hydrations and temperatures. The residence time for 
each region is calculated as 



= EE awE' 



(i) 



= 1 k=l 



where AUk is the fraction of time that the molecule i 
spent during the k-th visit in the region. For all hy- 
drations oscillations are present which are particularly 
marked for Nw = 500. The residence time reaches the 
maximum values at the surface. The structure of the resi- 
dence time profiles around the two main peaks is reminis- 
cent of the double layer structure of the density profiles 
(see Fig. |]). Apart for the layer attached to the Vycor 
surface the molecules have longer residence time in the 
regions where the density profiles have local maximal val- 
ues, while the residence time minimum is found roughly 
at the minimum of the density distribution between the 
two layers. The water molecules stay for a shorter time 
in the range between the two layers where the maximum 
of the water-water HB profile is located. So in the region 
of minimum density the HB formed are not very stable. 

The maxima of the residence time are enhanced by de- 
creasing the hydration. As discussed in Sec. IV this effect 
can be explained in terms of an increase of the rigidity 
of the hydrogen bond network at decreasing hydration. 

Close to the substrate where an elevated number of 
water- Vycor hydrogen bonds N^^ are present a com- 
parison of the average residence time in the zone 16 < 
R < 20 A with N^b shows that the two quantities are 
proportional in the range of temperature explored 



t res (T)=C N N™(T) 



(2) 



where the constant Cn increases at decreasing hydration 
level going from 10 for N w = 1500 to 20 for N w = 500. 

The distributions of the residence times in the different 
layers show peculiar features, as illustrated in Fig. [l2| 
for Nw = 1500 at room temperature. The distributions 
V{t rea ) for the layers more close to the surface are well 
described by a power law 



V(t) = at-" 



(3) 
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FIG. 10: Trajectory of the center of mass of a molecule in 
the xy plane d xy (t) = yj x 2 (t) + y 2 (i) for Nw = 500 and 
temperature T = 300 K (black line), T = 220 K (dashed 
line) . 



where for the molecules in layer 16 < R < 20 A /i = 1.51, 
while by considering the layer 14 < R < 20 A /i = 1.36. 
This behavior is a signature of a non brownian motion of 
the particles and it is observed at all the hydration levels 
explored in this work. Non brownian diffusion has been 
found recently in neutron scattering experiment and Mil 
simulation on water at contact with a globular proteinEll 
where the distribution of the residence times of the water 
molecules in the layers close to the protein surface is fitted 
to Eq. (|^) with exponents fi in fairly good agreement with 
the exponents found in the present work. 

For_water molecules far from the protein surface it is 
foundtj that the brownian regime is established with an 
exponential decays of the residence time distribution. In 
the inset of Fig. |l2| the residence time distribution is re- 
ported for the internal layer < R < 16 A of our system, 
also in our case the decay at long time can be fitted with 
an exponential 



V(t) = be 



-Bl 



(4) 



VI. SUMMARY AND CONCLUSIONS 



The structural properties of water confined in a cavity 
representing the average properties of the pores of Vycor 
glass have been studied by computer simulation. Water 
molecules are attracted by the hydrophilic surface and 
form a double layer structure of 4 4- 5 A. In this region 
it is found that the local properties of water are strongly 
dependent on the hydration level and the distance from 
the substrate. 

The hydrogen bond network is deformed with re- 
spectr±0|-bulk water, in agreement with experimental re- 
sults.LHaa At a distance greater than 4 A from the surface 
for the highest hydrations investigated (Nw = 1500 and 
Nw = 1000) the water molecules are arranged in a local 
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FIG. 11: Average residence time of the water molecules along 
the pore radius for T = 300 K (solid line) and T = 220 K (long 
dashed) at different hydrations: Nw = 500 (top), Nw = 1000 
(middle), N w = 1500 (bottom). 



tetrahedral order similar to bulk water, while approach- 
ing the substrate the water- Vycor hydrogen bonds sub- 
stitute the water-water hydrogen bonds and the tetra- 
hedral order is completely lost. At the lowest studied 
hydration (Nw = 500) the local order far from the sur- 
face is compatible with nucleation of a cubic crystalline 
phase at low temperature. Experimental evidences for 
nucleation of cubiCpice has been indeed found for water 
under confinement .Ej 

At a given hydration the average number of water- 
water HB increases with decreasing temperature while 
the number of nearest neighbors remains constant, this 
implies that the hydrogen bond network becomes more 
rigid upon supercooling. The effect is accentuated by 
decreasing the hydration. This findings are consistent 
with the recently proposed-Jragile to strong transition 
for strongly confined water .113 

The analysis of the residence time of the water 
molecules along the pore radius evidenced connections 
between this quantity and the density profiles, the wa- 
ter molecules persist for longer time in the regions of the 
maxima of the double layer structure and very close to 
the Vycor surface. 

The residence time distributions show an anomalous 
non brownian behavior in the layers close to the sub- 
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strate with a power law decay for long time. On the con- 
trary the distributions decay exponentially in the layers 
far from the surface. Experimental and computer simula=. 
tion studies on water at contact with a globular protcinEfl 
show a similar behavior. This agreement indicates that 
there are useful analogies between the behavior of water 
confined in Vycor and the behavior of water at contact 
with proteins. This is confirmed by a more refined anal- 
ysis of the dynamical properties of confined water in the 
low hydration regime.Ej 



FIG. 12: Residence time distribution at room temperature 
for N w = 1500 calculated for different layers: 16 < R < 20 A 
(triangles) and bottom curve, 14 < R < 16 A (circles) and 
top curve, in both cases the solid line is the power law fit 
(see Eq. ph. In the inset is reported the result for the layer 
< R < 16 A the solid line is the fit to the exponential law 
of Eq. (§). 
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